A g-deformed NJL model at finite temperature: 
chiral symmetry restoration and pion properties 

V. S. Timoteo 

Centra Superior de Educagao Tecnologica, Universidade Estadual de Campinas, 13484-370, Limeira, SP, Brazil 



oo 



C. L. Lima 

Instituto de Fisica, Universidade de Sao Paulo, CP 66318, 05315-970, Sao Paulo, SP, Brazil 

We review the implementation of a g-deformed fermionic algebra in the Nambu-Jona-Lasinio 
model (NJL). Lhe gap equations obtained from a deformed condensate as well as from the defor- 
mation of the NJL Hamiltonian are discussed. The effect of both temperature and deformation in 
the chiral symmetry restoration process as well as in the pion properties is studied. 
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O . I. INTRODUCTION 

o ■ 
(N ■ 

Chiral symmetry breaking is an important phenomena in hadron physics and is of fundamental importance for 
hadron properties. The difficulties involved in obtaining low-energy properties directly from QCD, the fundamental 
theory of strong interactions, have motivated the construction of effective models. Due to its simplicity and effective- 
ness in describing hadrons at low energies, the NJL model [l| has become the most familiar effective model for strong 
interactions. 

The standard scenario for spontaneous chiral symmetry breaking is the arising of a quark condensate, fn the NJL 
I , model the condensates appear when the strength of the contact interaction exceeds a critical value, separating the 
■ Wigner- Weyl and Nambu-Goldstone realizations of chiral symmetry. The condensate < qq > is the order parameter. 
Q\ \ Apart from giving a reasonable description of the light mesons, the NJL model is also very useful for studying the 
chiral symmetry breaking process as well as its restoration at finite temperature 0, HI • The NJL model is in addition 
very suitable for testing new ideas. In previous works, we have investigated the effect of g-deformation in the chiral 
symmetry breaking process within the context of the NJL model. In particular, we observed how the condensate is 
affected by the deformation. As a direct consequence, the dynamical quark mass, the chiral symmetry breaking and 
restoration processes are accordingly affected 0, 0| . 

The aim of this work is to study the effect of both temperature and deformation in the chiral symmetry restoration, 
and evaluate the pion properties in a g-deformed NJL model at finite temperature. This paper is organized as follows. 
1 ■ Sec. II introduces the g-deformed fermionic algebra which will be used along this work, and shows how the deformed 
algebra is implemented in the NJL model. In Sec. Ill, we present the investigation of the chiral symmetry restoration 
^ , with both finite temperature and g-deformation. Finally, Sec. IV contains our main conclusions. 
Ch ■ 
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• r-j . II. g-DEFORMATION IN THE NJL MODEL 

X' 

The Nambu-Jona-Lasinio model was first introduced to describe the nucleon-nucleon interaction via a four-fermion 
contact interaction. Later, the model was extended to quark degrees of freedom becoming an effective model for 
quantum chromodynamics. Details of the NJL model and the Bogoliubov-Valatin variational approach can be found 
in review works 0, 13 ■ 

In this section we discuss the implementation of the deformed algebra to the NJL model. The g-deformed fermionic 
algebra Q that we shall use is based in the work of Ubriaco where the thermodynamic properties of a many 
fermion system were studied. In the construction of a g-covariant form of the BCS approximation 9], it was shown 
that the creation and annihilation operators of the su q (2j + I) fermionic algebra are given by 

A im3 = a jmj Yl ( 1 + Q a )i a 3<) » (!) 



<3\ 

o 
in 
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Am, = 4m t II i 1 + Q a lflo) » ( 2 ) 
i—rrij +1 

where Q = g _1 — 1, j = 1/2 and rnj = ±1/2. As discussed in the first consequence of the above deformation is 
that only the operators corresponding to rrij — — ^ are modified, meaning that only negative helicity quarks (anti- 
quarks) operators will be deformed since in the NJL model we deal with quarks (anti-quarks) creation and annihilation 
operators. 
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We have two different approaches to obtain a new gap equation. The first one consists in to perform the q- 
deformation in the condensate, which is a part of the self-consistent equation for the dynamical mass. In the second 
one the g-deformation is performed in the NJL Hamiltonian, and we use the Bogoliubov-Valatin procedure to obtain 
the new gap equation. We now turn to a detailed discussion of these two different approaches. 



A. Deforming the condensates 

To obtain the deformed gap equation we work with the BCS-like vacuum in the standard Bogoliubov-Valatin 
variational approach 

\NJL)= Yl [cos9{p) + ssm9(p)tf(p,s)J(-p,s)] |0) , (3) 

p,s=±l 

which, for a given momentum p, is expanded as 

\NJL) = cos 2 9{p) |0) 

+ sin 9{p) cos 0(p)&t ( p , +)d t (-p, +) |0) 
— sin 9{p) cosO(p)tf(p, -)dt(-p, -) |0) 

-sin 2 0( J5 )6t(p,-)rft(_p ; _) 6 t(p, + ) o jt ( _p i+) |o) . (4) 



The quark fields are expressed in terms of g-deformed creation and annihilation operators as 

„ J \2tt) 



(5) 



where the g-deformed quark and anti-quark creation and annihilation operators B, B\ T>, and T>\ are expressed in 
terms of the non-deformed ones according to Eqs. (JIJ and J2J 

B- = 6- (l + Qb\b^j , Si = &i (l + Qb + b+^j , (6) 

2?- = d_ (l + » p - = d - ( x + Q d + d +) > ( 7 ) 

£?+ = 6+ , Bl = b\, (8) 
2?+=d+ , 23+. = (9) 

where + (— ) stands for the positive (negative) helicity and the notation has been simplified: (p, s) — > s. We would like 
to note that, as discussed in Ref. (3, the deformed vacuum differs from the non-deformed one only by a phase and, 
therefore, the effects of the deformation comes solely from the modified field operators. Additionally, the q-deformed 
NJL Lagrangian, constructed using ip q instead of ip, is invariant under the quantum group SU q (2) transformations. 
This can be seen by using the two-dimensional representation of the SU q {2) unitary transformation given in Ref. 0. 
The deformed gap equation is 

m = -2G(jH>) q , (10) 

where (tpip) is the g-deformed condensate calculated using the BCS-like vacuum Eq. 0, 

(^) q = (NJL\^ q \NJL) 

= (W) + {NJL\Q\NJL), (11) 

where (ipip) is the non-deformed condensate and {NJL \ Q\ NJL) represents all non- vanishing matrix elements arising 
from the g-deformation of the quark fields. The contribution of these g-deformed matrix elements is 

{NJL \Q\ NJL) = Q f [sin 26{p) - sin 29{p) cos 29{p)\ . (12) 

J (2n) 
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The calculation of the deformed condensate will be performed in a similar way as in the usual case. It requires also 
a regularization procedure since the NJL interaction is not perturbatively renormalizablc. For reasons of simplicity a 
non-covariant trimomentum cutoff is applied arising 



<^> 9 



3to 



dp 



P 



Q 



Vp 2 + ™ 2 2 



dp- 



P 



(13) 



for each quark flavor. At this point we see that the dynamical mass is again given by a self-consistent equation since 
the condensate depends also on the mass. Inserting Eq. (jT3")l into Eq. (fTTTfl we obtain the deformed NJL gap equation 



2Gm 



m 



1 



dp 



Q 



Vp 2 + < 



dp- 



It is easy to see that for Q = (q = 1), we recover the NJL gap equation in its more familiar form 

,2 



P 



2Gm f A , 
m = — — / dp- — 
/o VP 2 + : 



7T^ 



to , 



(14) 



(15) 



where too appears only if we consider the current quark mass term C ma ss = —mot/ji/j in the NJL Lagrangian. 

The pion decay constant is calculated from the vacuum to one pion axial vector current matrix element, which, in 
the simple 3D non-covariant cutoff we are using 0, is given by 



fl = N c m 2 



d 3 p 



1 



^3/2 > 



(16) 



for each quark flavor. The deformed calculation of /„■ is performed directly by substituting the dynamical mass in 
Eq. (|TB1) from the one obtained in Eq. ijLT) . 

As in the non-deformed case, the (/-gap equation has non-trivial solutions when the coupling G exceeds a critical 
value G c related to the cutoff. The dynamical mass is accordingly modified through the deformed gap equation i|14fl . 

The behavior of the condensate around the critical coupling, G c , is similar for both deformed and non-deformed 
cases, meaning that the adopted procedure to g-deform the underlying su(2) algebra in a two flavor NJL model does 
not change the behavior of the phase transition around G c . The formalism developed in |8j,|9j] allow q- values smaller 
than one (which corresponds to Q > 0). It is worth to mention that in this case the g-deformation effect goes in the 
opposite direction, namely, the condensate value and the dynamical mass decrease for q < 1. 



B. Deforming the NJL Hamiltonian 



Here it is important to stress again that the deformed version of this vacuum differs from the non-deformed one 
only by a phase and, therefore, the effects of the deformation comes solely from the deformation of the Hamiltonian. 

The deformed functional for the total energy will be obtained from the vacuum expectation value of the q-deformed 
NJL Hamiltonian: 



W [6{p)] - (NJL\H q NJL \NJL), 

where the g-deformed Hamiltonian is 



(17) 



(18) 



and ip q is given by Eq. JSJ). Due to the additive structure of the g-deformation, the deformed Hamiltonian can be 
written as 



Wnjl — Hnjl + W(Q), 



(19) 



and the functional will have the same form 



W q [0(p)} = W[6(p)}+W[Q,6(p)] 



(20) 
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The last terms of Eqs. (|19|) and Q20J1. namely H(Q) and W[Q,0(p)], stand for the new terms of first order in Q 
generated when the algebra is deformed, and, therefore, they vanish for q = 1 (Q = 0). 

The non-deformed case, which corresponds to q = 1 (Q = 0), has been discussed in Sec. II. In order to obtain 
the new gap equation, we need to calculate the new matrix elements arising from the g-deformation of the NJL 
Hamiltonian, and add them to the non-deformed functional. This procedure yields to the full g-deformed functional 
for the total energy: 



d 3 P ± , Ar Ar Q f d 3 p 2 

2 



W q [6(p)} = -2N c N f {l + Q) / — JL P cos <l> + N c N f f / — ^cos : 



(27r) J 1 2 J (2tt) 

+N c N f ^QGA 3 J ^co S <f>-(N c N f fG(4 + Q)U -^smcfi} . (21) 



Defining the new variables 



P q = P + P , (22) 
P = (1 + Q)p, (23) 

6tt 2 



p o = —^fGK 3 Q , (24) 



* = -f?. ( 25 ) 
G = g(i + Q) , (26) 



4, 

and performing the same minimization procedure as in the non-deformed case we obtain 

P q tan (f> + K sin <j> = 4G'N c Nf J ^-y-sin^. (27) 

The term Ksvcicf) was overlooked in a previous paper p| [l^. Trying to keep as much similarity as possible with 
the usual solution, instead of solving the equation P q tan(j)(p) + Ksm<f>(p) = M, the sin cf>(p) will be substituted by 
atan0(p), with a = 0.59 chosen to minimize the difference between sin^> and tan</> in the interval [0, 7r/4]. 
The new gap equation then becomes 

M = AG' N c Nf [ , = , (28) 

f J (2^) 3 y^ + M 2 

provided the variational angles have the same old structure (but they now are g-dependent) 



tan 29 q (p) = , sin 20, (p) = , cos 20, (p) = ' g , (29) 

^« ,/K 2 + M 2 ,/K 2 + M 2 



M M K n 

>q\v) - — ,sm20,( 

K 

Vg | ^ y 1 

where 

K q = P q + aK. (30) 

It is easy to see that, when q — > 1 (Q — > 0), Eqs. (|27|l . (|28|l . and i|29|) reduce to their non-deformed versions, since 
K q -> p and G' -> G. 

In analogy with the non-deformed case we can write the gap equation in terms of the quark condensates as 

M = -2G . (31) 

Comparing the two forms of the gap equation Eqs. ((14(1 and 1)31(1 . we find a new deformed condensate given by 

K 2 + M 2 
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This condensate is different from the one obtained in the previous section, where the condensate was explicitly 
deformed [4|. It also has exactly the same form of the non-deformed one, but is written in terms of the new variables. 
It is worth to mention that the new condensate is not obtained by calculating the vacuum expectation value of a 
deformed scalar density. In fact, it corresponds to the gap equation which arises from the variational procedure started 
from the g-deformed Hamiltonian. The new pion decay constant can also be obtained in analogy to the non-deformed 
case 

Ft = NcM*[ A ^L 1 (33 ) 



W ( K 2 + M*) 



A closer look to the left side of Eq. shows that the dynamical mass has two components: one proportional to 
the original dynamical mass m and another term which we call M 

M = [1 + (1 - a/2)Q] p tan 26 (p) + p tan 26{p) , (34) 

rn Mn 
v » ' 

M dyn 

where po — — ^ GA is Pq (Eq. 124ft for N c = 3 and Nj — 2. We can then calculate the effect of the deformation on the 
dynamical mass m: 

m = l + (l-a/2)Q ' (35) 

where M is obtained by solving Eq. 131|) . In an analogous way, the effect on the condensate can be calculated by 
substituting M — ► m in the numerator of Eq. Q32J1 t l3| • 



III. TEMPERATURE AND q-DEFORMATION 



In this section we want to study the interplay of temperature and q-deformation in the NJL model. We review the 
standard approach to introduce finite temperature and chemical potential in the NJL model 0,0 in subsections IV. A 
and IV. B. When temperature is introduced in the model, the condensate is replaced by the thermal expectation value 
of the scalar density, which contains the Fermi-Dirac distributions. We also have a gap equation for the chemical 
potential, which is modified by the interactions in such a way that we need to solve a system of coupled gap equations. 
In subsections IV. C and IV. D, the standard formalism is extended to incorporate the effects of q-deformation. In 
particular, aspects of chiral symmetry restoration and pionic properties in the q < 1 regime will be discussed. 



A. The NJL model at Finite Temperature 



The starting point for the study of the thermodynamics of the NJL model is the partition function 

Z = Tr exp [~(3 (U NJL - MiM)] , (36) 

where = T _1 , Mi is the valence quark number operator of flavor i, and /ii are the corresponding chemical potential. 
From the above partition function, one can calculate the thermal expectation value of an operator A 

{(A)) = ^TrA{exp[-/3 (H NJL - p l M t )}} , (37) 

where the operator in question can be, for example, the scalar density tptfj or the quarks density ip^ip. 
Here the Hamiltonian is given by 

Hnjl = -iip-f ■ V?A - G (ipip) 2 - G ftjk'ysTip) 2 + m ip^, (38) 
= -i-ipj ■ Vip - L int + m ip->p, (39) 

which, in the mean field approximation, is written as 

Q Q 

Hmf = -vh ■ + mtpip + Gal + —cfzM - jjj^of, ( 40 ) 
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where m is the effective quark mass 



M is the number operator 



and i7i, C2 are defined as 



m = m,Q-2Gax, (41) 

TV = ^70^ = tp^tp, (42) 

*i = ((^» , (43) 

<T2 = ((^7o^»- (44) 

The Hamiltonian in the mean field approximation represents a system of free particles of mass m and the chemical 
potential is given by 

G 

M = Mo--^-0"2, (45) 

where fj,o is the chemical potential when there is no vector interaction, generated by the Fierz transformation of Ci n t- 
Equations l|41|l and (|45|l form a system of self-consistent coupled equations 

(m = m -2Ga 1 , 
\ fi = - ^cr 2 . 

Solving the above equations we obtain the effective mass and chemical potential, which we can put back in the 
expressions for <j\ and <72 to obtain the condensate and the density at a given temperature T. 

B. Thermal Expectation Values 

The dynamical mass calculation depends on expectation values containing the Fermi-Dirac distributions. The 
thermal expectation values relevant for the NJL model, namely a\ and 02, can be calculated within the formalism 
of the thermal Green function [TlT | , which has the advantage of separating the parts that depend on the temperature 
and the chemical potential. The thermal Green function for a free fermion at a temperature T and chemical potential 
fi is written as 

S F (p; T, fi) = 2 * + ™. + 2mS (p 2 - m 2 ) 

x (j> + m) [6 (p°) n(p, T,fi) + 8 (-p°) n(p, T, M )] , (47) 



where 



and 



1 + exp [—p (E - n)\ 



"(P. T ' M) = T- F fl , F , yj , (49) 

1 + exp [—/? (Js + /i)J 



are the fermions and anti-fermions distribution function respectively with E = -\/p 2 + m 2 . 

Making use of the anti-commutation relation for the fermionic fields, the expectation values can be written as 

((^(x)V(x)» = -Tr lim «T [^(a#(y)]>) , (50) 

((? 7 o^» = -Tr hm 7o «T ty(a#(y)]» . (51) 

The term on the right-hand side of the above equations is the definition of the thermal Green function in the 
configuration space: 

iS F (x - y- T, (i) = ((T )) . (52) 
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The expectation values are then re-written as 



= -i lim Tr f -^S F (p;T, e^*"^, (53) 

_/ (27r) 

«?7o^>)) = -» Hm Tr / -^7oS F (p; T, A*) e^*"^. (54) 



Applying a cutoff in the momentum, the results for a\ and (72 are the following 

N c N f f A 2 m 



a, = = V / d PP 2 ^ [1-n-n, (55) 

^ Jo & 

°2 = «?7o^» = / <fcp 2 [n - n] , (56) 



7T- 



(I 



and the barion density can be written as 

PB = \ ((^7oV>» = ^ jf ^ ["-»]■ ( 57 ) 

C. Condensates, mass, and chiral symmetry restoration 

An interesting window to the application of g-deformation to hadronic physics can be opened for g-values smaller 
than one. It is tempting to explore the behavior of the condensate in this new regime for smaller values of q, even 
considering that the truncation at order Q may not be granted. In this case, we observed in Q that the chiral 
symmetry is restored in the limit q — > 0, since the condensate vanishes. It is then worth to investigate the effect of 
both temperature and g-deformation in the chiral symmetry restoration process as well as in the pion properties. 

In order to study the effect of both temperature and g-deformation in the chiral symmetry restoration process, we 
can replace the gap equation for the mass in the system of coupled equations by its deformed version Eq. (|58|) 

M = m - 2G' ((##)) , (58) 

so that the new system of coupled equation we need to solve is 

M = mo-2G'53i, 

G (59) 

P = PO- ]^2, 

where Si represents the deformed version of the thermal condensate <j\. 

Solving the new set of coupled gap equations we can observe the condensate as a function of temperature for 
different values of the deformation parameter. It is important to mention that the effects of the algebra deformation 
come from gap equation for the mass. However, the chemical potential is also affected once it depends on the mass 
through the coupled equations. Thus, for a given temperature, we obtain the mass and the chemical potential by 
solving the system of coupled gap equations. We then calculate the condensate and the density. 

In order to visualize the chiral symmetry restoration process with the influence from both temperature and defor- 
mation, we plot the surfaces shown in Fig. ^where- the chiral condensate and the constituent quark mass are plotted 
as functions of both temperature and g-deformation. 

At q < 1 deformations, an interesting result is observed: the chiral symmetry restoration is similar to the one 
observed by Klimt et al. as function of temperature, T, and density, p. However, the process is slower as q 
gets smaller when compared with the case of p getting larger. Figure H] shows that the condensate is substantially 
affected by the deformation and lowers as q decreases, while the effect of temperature is softened. The same behavior 
is observed in the dynamical quark mass since it is linearly constrained to the condensate by the gap equation. 

We have also calculated the barion density as a function of T and q, which is displayed in Fig. [21 The effect of 
the deformation on the barion density in only observed as the temperature increases. The density also lowers as q 
decreases, and the temperature effect is also smoothed down. 

The behavior of the dynamical mass and the chiral condensate directs to chiral symmetry restoration. We would 
therefore expect the pion to become massless in the limit q — > 0. In the next section, we evaluate the pion decay 
constant and mass in order to investigate the pion properties in a regime where condensate is lower than in the 
standard NJL model due to the effects of the g-deformation. 
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D. Pion Properties in the q < 1 regime 



Recent works |l3j have shown that the spontaneous chiral symmetry breaking may result from a balance of the 
density of states and their mobility in the medium. It may occur with a vanishing < qq > condensate as long as the 
low density of states is compensated by a high mobility. In this case, the condensate is no longer the order parameter. 

The effective models for QCD, like the NJL, fail to describe hadronic observables when the condensate is small or 
null. The g-deformed version of the NJL model seems to be an alternative to study the low condensation scenario. 
As a first application, we start with the pion mass since it cannot be reproduced in the standard NJL model in a low 
condensation regime. 

With a fixed value for the current quark mass, after obtaining the pion decay constant through Eq. 1331 we obtain 
the pion mass by using the Gell-Mann-Oakes-Renner (GOR) relation 



written in terms of our deformed quantities. The pion properties (decay constant and mass) as functions of both 
temperature and q-deformation are shown in Figure [3J In order to have an idea of the effect of the deformation in a 
low condensation regime, we consider the case of q = 0.5, where condensation is approximately 35% lower than in the 
non-deformed case (see Fig. QJ. In this particular scenario, we obtain a pion mass of about 190 MeV. The standard 
NJL model would give a pion mass of about 100 MeV in such condensation scenario. The correlations between the 
constituents of the system introduced when the underlying su(2) algebra is deformed seems to compensate the low 
condensation, which is responsible for keeping a reasonable mass for the pion even when the value of the condensate 
is getting smaller. Recent measurements of the semileptonic decay K + — » ir + ir~ e + v e ^3] obtained a value of the 
hit S-wave scattering length which indicates that the condensate is indeed the order parameter |15| . However, it is 
interesting to note that the deformed version of the NJL model is able to describe the pion even with smaller values 
of the condensate. 



So far we have performed the ^-deformation of the NJL model with finite temperature also taken into account. We 
studied the g-deformation of the underlying su(2) algebra in a two flavor version of the NJL model and investigated 
an important feature of chiral symmetry: its restoration at finite temperature. Also, we studied the pion properties 
in a low condensation scenario, where the standard NJL model would underestimate the pion mass. 

As far as the effects of both temperature and deformation are concerned, our main conclusions can be summarized 
as follows. The quark condensate and the dynamical mass decrease as q gets smaller than one, and chiral symmetry 
would be restored in the limit q — » This effect is specially observed at low temperatures, since we have thermal 
restoration at large T. Surprisingly, as the condensate lowers due to the deformation, the pion not only keeps massive 
but also increases. The GOR relation, Eq. 1601 with our deformed quantities gives the reason for this behavior, where 
the pion decay constant decreases faster than the square root of the condensate. 

It seems that the effect of the deformation compensates the lower condensation, suggesting that the quarks' 
mobility is affected when the algebra is deformed. In fact, the correlations introduced by the deformation generate an 
extra momentum Pg for the constituents of the system (see Eq. I24[) which is responsible for enhancing the mobility 
in the medium. 
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FIG. 1: The chiral condensate and the dynamical quark mass as functions of both temperature and g-deformation in the q < 1 
regime. The results presented in the figures were obtained with the following parameters: G = 6.58 GeV~ 2 , A = 0.6 GeV, 
m = 5.87 MeV, and ^ = 0.35 GeV. 
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FIG. 3: The pion decay constant and the pion mass as functions of both temperature and g-deformation in the q < 1 regime. 



